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ABSTRACT 


We present in this study a three step algorithm for the decomposition of 
arbitrary, three-dimensional, planar polygons into convex polygons. Through a 
series of translations and rotations. an arbitrary polygon is mapped onto the x-y 
plane. then broken into a set of convex polygons, and finally mapped back to the 
polygons original coordinate system for filling and display by special graphics 


hardware. 
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I. INTRODUCTION 


Computer hardware advancements applied to the field of computer graphics 
have recently made possible the real-time display and animation of three- 
dimensional object models. In most currently produced graphics systems, the 
three-dimensional (3D) object models are comprised of sets of 3D polygons. For 
vears the key bottleneck to generating real-time images of 3D models has been the 
limited speed with which the 3D polygons comprising the model have been filled. 
This bottleneck has been alleviated recently with the advent of special hardware 
designed to fill polygons at rates exceeding 44 million pixels per second [Ref. 1]. 
For real-time animation of 3D object models, such polygon fill rates are essential. 

The majority of the graphics systems that provide polygon fill hardware only 
provide hardware for filling convex polygons. Concave polygons passed to such 
hardware usually are filled incorrectly, if at all. To utilize the capabilities of such 
hardware for the general polygon case, that is, a mix of concave and convex 
polygons, it is necessary to have some means of decomposing an arbitrary 3D 
planar polygon into a set of 3D planar convex polygons completely covering the 
area defined by the original polygon. The following study attempts to provide 
such an algorithm. We then discuss the implementation and performance 


characteristics of that algorithm on the Silicon Graphics, Inc. IRIS workstation. 


Il. BACKGROUND 


A. DEFINITIONS 

A polygon is defined as a geometric figure that is described by an array of x, 
y. and z coordinates for its vertices. A two dimensional polygon lying in the x-y 
plane is a special case of a three dimensional polygon in which the z coordinates 
are all zero. A convex polygon is one that has no points of negative curvature 
along its boundary [Ref. 2: p. 217]. This means that a line can be drawn from 
any point inside the polygon to any of its vertices and have the line lie entirely 
within the polygon boundary [Ref. 3: p. 350]. A concave polygon fails the above 
test. i.e. we cannot draw a line from every point inside the polygon to every 
vertex and have the line remain within the polygon boundary. However. any 
point inside a concave polygon can be shown to lie inside a convex polygon 
created by some of the vertices and sides of the original polygon [Ref. 3: p. 330]. 


Figure 2.1 provides examples of convex and concave polygons. 


B. DECOMPOSITION 

To fill a concave polygon utilizing hardware that only correctly fills convex 
polygons, we must first break down the polygon into simpler convex polygons that 
together cover the area of the original. The process of decomposing a concave 


polygon does not necessarily yield a single unique set of convex polygons, but the 


Nod © 


convex polygons 


AS @ 


concave polygons 





Figure 2.1 - Convex and Concave Polygons 
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Figure 2.2 — Two Equivalent Decompositions 
of a Concave Polygon. 





result of two different decompositions of the same polygon should be equivalent as 
shown in Figure 2.2. 

During the operation of decomposing a polygon. there are two main checks 
that must be conducted to ensure complete and proper decomposition takes place. 
The first check is to ensure that no area of a created convex polygon falls outside 
the boundary of the original concave polygon. The second check is to ensure that 
the concave polygon is completely covered by the total areas of the created convex 
polygons. Figure 2.3 provides examples of the two checks. 

The polygon decomposition algorithm produced for this study is comprised of 
three basic steps. The first maps the three dimensional coordinates of the polygon 
into a two dimensional plane. This process generates a polygon that is more 
easily examined and decomposed. The second step conducts the actual 
decomposition an performs the two decomposition checks. The third step maps 
the coordinates of the created convex polygons back to the the original orientation 
in three-space. At this point. the created convex polygons are then ready for 


passing to the graphics system’s polygon fill hardware. 


C. RESEARCH FACILITIES 

The IRIS (Integrated Raster Imaging System) 2400 Graphics Workstation, 
manufactured by Silicon Graphics, Inc., is the target system for our polygon 
decomposition algorithm. By incorporating custom VLSI chips to replace slower 


graphics software. the IRIS system can accept polygons in_ user-defined 
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Figure 2.3a - Without a check to ensure the convex 
polygon does not fall outside the bounds of the 


original polygon, vertices could be chosen that 


yieid an incorrect decomposition. 
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Figure 2.3b - The second decomposition check ensures 
the concave polygon is completely covered by the 


decomposed polygons. In this case all vertices 
and edges were utilized, however a hole is still 
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coordinates, and transform them to screen coordinates for display on a color 
graphics monitor that provides 1024 x 768 resolution. This use of special 
hardware yields processing speeds that are some 100 times faster than similar 


operations carried out in software. 
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III. STEP I: MAPPING THE 3D POLYGON 0 THE@D ey Pi 


A. THE NORMAL 1O 4A. POerveo. 

The first step of the polygon decomposition algorithm, is comprised of several 
small steps. The goal of these steps is to map the polygon to the x-y plane. First. 
three non-colinear points are selected from the array of x, y. and z coordinates of 
the polygon. These are used to define the equation of the plane for the 3D 


polygon. and the normal to that plane (Figure 3.1). 


Ax + By +.¢€z —) 
Ne A. + B; + Cy 
The coefficients A. B. and C are calculated from the coordinates of the three 
selected points by the following equations: 
A = (Yo — Yi) (Z2 — 21) — (¥2 — Y1) (Zo — 21); 
B = (Zp — 21) (X2 — &) — (22 — 701 X,): 


C = (Xp — X1)(¥2 — yi) — (X2 — X1) (Yo — Y1)- 


The coefficients are then used to solve for D of the equation of the plane. 


12 





Three vertices O, 1, and 4 
are selected to determine the 
Normal to the plane. 


Figure 3.1 The Normal to a Polygon 
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B. ANGLES ABOUT THE AXES 

Once the equation of the plane is determined. each point in the polygon input 
array is checked to ensure that the polygon is in fact planar. If the polygon 
passes this test, then the angles between the normal and the axes of the 
coordinate system are found. From these angles, the necessary rotations. about 
the x and y axes, to bring the polygon into the x-y plane, are calculated. The 
following equations [Ref. 4: pp. 24-59] are used to calculate the angles shown in 


Figure 3.2a. 








eee 
\/ Agee °C 
hag 
J A* + B*+C 
6 
cos 7 = 





JA*+B*+C 


For the polygon to lie in the x-y plane, the angle between the normal and the 
Z axis is either zero or 180 degrees. Necessarily then rotations are conducted to 
bring the angles between the normal and the x and y axes to 90 degrees each. 

First a projection of the normal onto the x-z plane is made as shown in Figure 
3.2b. The dot product of that projected normal \’ and the z axis is calculated. 
The dot product is divided by the product of the magnitudes of the two vectors 


yielding the cosine of the angle of rotation about the y axis (4). 
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Figure 3.2 Angles about the Axes 
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A second projection of the normal onto the y-z plane is made and the angle of 


rotation about the x axis (@) is calculated in the same manner as above (Figure 


Sec): 


C. TRANSFORMATION MATRIX 

To actually rotate the polygon into the x-y plane, a transformation matrix 
must be created that is used in the multiplication of the polygon’s coordinates. 
We first translate the polygon to the origin by adding the negative x, v, and z 
coordinates of one of the polygon’s vertices. Assuming 1, m. and n are the 


coordinates selected. the matrix becomes: 


1 O OO 0 
0 1 0 O 
0 O 1 O 
—~] -m -n 1 


A rotation matrix for the angle ¢ about the y axis is then created 


cos@ O sing O 
f) 1 O O 
—sing O cos¢ 0 
0 oO oO 1 


and concatenated with the translation matrix. This is followed by a multiplication 


with a rotation matrix for the angle @ about the x axis. 
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1 O 0 O 
O cos#@ —sin@ O 
O sin8@ cos@ O 
0 QO 0 1 


The result of these matrix concatenations is the following transformation matrix 


that is used to map the array of vertices onto the x-y plane. 


COS —sin @sin8 —singcosé 0 

0 cos8 sind 0 

sing —cos¢sin8 cos¢cosé 0 
—lcos¢@—ncos¢ —mcos6+ (Ising—ncos¢)(—sin#@) —msin#@+(Ising—ncos¢)(cos8@) 1 


Once we have this matrix set up. call it M,,ine, We then use it to transform 
the polvgon’s coordinates via a simple row vector by matrix multiplication. The 
result of this multiplication is a transformed set of coordinates. in which the z 
coordinates are zero. The x and y values of this transformed array are then used 
in the decomposition step of our algorithm. We note at this point that for the 
mapping of the coordinates back to the original vertices (Step III of the Alg.), 
that the index into the transformed array is also the index for the untransformed 


point in the original array of vertices. 


1h 


IV. STEP Il: POLYGON DECOMPOsSTiliG x 


The decomposition step divides the polygon into multiple convex polygons. 
In so doing. the algorithm performs the two checks to ensure that no area outside 
the polygon is shaded and that the created convex polygons completely cover the 
area of the concave polygon. This step requires that the vertices of the polygon 
be ordered in a clockwise fashion. with respect to the +z directed normal to the 
x-y plane. Thus a determination of the ordering must be made and if necessary 


the array of vertices reversed. 


A. CLOCKWISE ORDERING 
To determine if the vertices of the polygon are ordered clockwise. the 
polygon’s area is calculated using the formula [Ref. 5: p. 369] below. (Note z,,y, 


are the polygon’s 2D coordinates. ) 


= ((xoy1) + (x12) + irs + (X,-1Yn) +(XnYo) —(YoX1) —(¥1X2) — 7 — (¥n-1%,)— (aaa 


The above equation yields a negative area for clockwise and a positive area for 
counterclockwise polygons. If the polygon is ordered counterclockwise then the 
array of x, y, and z coordinates that define the polygon is reversed before the 


selection procedure is begun. 
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B. VERTEX EXAMINATION 

The selection of vertices for a convex polygon is conducted by examining each 
vertex of the concave polygon with respect to the previously selected vertices of 
the developing convex polygon. This is begun by initially selecting two adjacent 
vertices. then in clockwise order, examining each of the polygon’s remaining 
vertices to build as large a convex polygon as will fit within the boundary of the 
concave polygon. 

A series of tests are conducted in the selection of vertices for the convex 
polygon. As vertices are selected. the line segments connecting these become the 
edges of the convex polygon and define its boundary. 

1. Test One 

The first test of a vertex is if it lies to the left or right of each line 
segment of the developing convex polygon. The vertices of the concave polygon 
are examined in clockwise order. Thus if a vertex lies to the left of any of the line 
segments of the developing polygon, then a concave angle is necessary to include 
the vertex in the polygon. As shown in Figure 4.la, vertex 2 lies to the left of the 
line segment 0-1. A concave angle ( » in Figure 4.1b) is necessary to add vertex 2 
to the developing polygon. Therefore vertex 2 is not selected and the process 
moves on to perform the same test on vertex 3. 

2. Test Two 

A vertex that passes the first test may still require a concave angle to be 


included in the developing convex polygon. In Figure 4.2. the vertices 0, 1, 2, 3, 
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Figure 4.1 - (a) Vertex 2 lies to the left of the 


line segment O-1, (b) Concave angle () would 


be required to connect Vertex 2 to the devel- 


oping polygon, (c) Vertex 2 is skipped and 


the examination is continued with vertex 3. 
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(f) 





Figure 4.2 -— (a) Vertex 5 lies to the right of each 
line segment of the developing polygon, (b) Concave 
angle (p) is required to connect Vertex 5 to Vertex 


O, (c) Concave angle (pf) can be identified since O 
lies to the left of line segment 4-5, (d) Vertex 4 


is skipped and the test repeated with line segment 
3-5, (e) Vertex 3 is skipped and the test repeated 
with line segment 2-5, (f) Result of the tests are 


that vertices O, 1, 2 and 5 are selected. 
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and 4 have been previously selected for a convex polygon with vertex 5 under 
examination. Vertex 5 passes the first test since it lies to the right. or inside of 
each of the line segments. However, a concave angle (p in Figure 4.2b) is again 
required to connect vertex 5 with vertex 0 of the developing polygon. A concave 
angle of this type is identified, as shown in Figure 4.2c. if the initial vertex of the 
developing convex polygon lies to the left of the line segment created bv the 
vertex under examination and the last selected vertex of the developing polygon. 
In Figure 4.2c. the initial vertex. 0. lies to the left of line segment 4-5. To remedy 
this. the last selected vertex of the developing convex polygon, vertex 4, is 
removed and the test is again conducted with the line segment formed by the 
vertex under examination and the last selected vertex (Figure 4.2d & 4.2e). This 
process is repeated. removing vertices from the developing polygon, until the 
initial vertex no longer lies to the left of a new line segment. 
3. Test Three 

The third test conducted is to ensure that no line segment created in the 
convex polygon intersects any of the edges of the concave polygon. Such an 
occurrence results in area outside the concave polygon being included in the 
convex polygon and an inaccurate polygon fill taking place. In Figure 4.3, the 


ad 


previously selected vertices are 0. 1. 2, and 3. Vertex 7 is under examination and 
passes the two earlier tests. However. the selection of vertex 7 results in the 


convex polygon 0, 1, 2, 3. 7 which contains area outside the concave polygon. A 


situation such as this is identified since line segment 3-7 intersects line segment 
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Figure 4.3 - Tests must be conducted for the 


intersection of developed line segments with 


the edges of the original concave polygon. 
Otherwise Vertex 7 would be an acceptable 


vertex for inclusion in the developing 
PommeonsOr 1° 2° 3S... 
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13-0 of the concave polygon. When this occurs. the vertex under examination is 
removed from the convex polygon and the next vertex of the concave polygon is 


examined. 


C, SELECTION COMPLETION 

The selection process for the developing convex polygon concludes when all 
vertices of the concave polygon have been examined. The convex polygon is then 
ready for the third step of our algorithm. mapping to the original polygon’s 3D 
coordinates. Figure 4.4 is an example of a concave polygon that was decomposed 
to the three shaded convex polygons using the above three tests. We have a 
"hole" left in the concave polygon even though all vertices and edges of the 
concave polygon have been utilized. Holes such as this are prevented by taking 
the initial and last vertices of a developed convex polygon as the starting vertices 
of a convex polygon to be developed. In Figure 4.4, vertices 0 and 2 of polygon 0, 
1, 2 are taken to develop convex polygon 0, 2, 4 to fill in the "hole". Once this is 
accomplished, step two of the algorithm is repeated by selecting a vertex of the 
concave polygon that has not yet been utilized in a convex polygon. The 
decomposition process is complete when there are no vertices or edges of the 


concave polygon that have not been utilized in one or more convex polygons. 


24 





Figure 4.4 - All vertices and edges of 


the polygon have been used, yet a hole 


remains. 
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V. STEP Ol: MAPPING TO THE ORIGINAL COORDINATE SYSTEM 


The second and third steps of the decomposition process are interleaved, i.e., 
each created convex polygon is immediately mapped back to the original 
coordinate system. This is accomplished by maintaining the original array of 
vertices that the user provided for the three-dimensional polygon. As a convex 
polygon is produced. it is actually represented by the indices of its vertices. Actual 
vertex coordinates are extracted from the vertex array. as required. using the 
index. In this manner. the transformed coordinates can be utilized during 
decomposition and the index can be used to obtain the original coordinates for 
final display. By this method, there is no consequential need to translate and 


rotate the vertices back to the original coordinate system. 
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VI. IRIS IMPLEMENTATION 


Implementation of the decomposition algorithm on the IRIS system is via 
approximately 1600 lines of C code (Appendix). Several data structures and 


various IRIS system hardware features are utilized as described below. 


A. MATRIX IMPLEMENTATION 

In the development of the transformation matrix, IRIS system calls for 
translation and rotation are used to alter the top of the system matrix stack. 
These routines perform the matrix concatenations needed to translate and rotate 
the polygon to the x-y plane. To convert the polygon vertices to 2D coordinates, 
four vertices at a time are loaded into a matrix. The x, y, and z coordinates go 
into the first three columns with 1’s filling the fourth column to provide the 
necessary homogeneous coordinates. This is placed on the system’s matrix stack 
and multiplied with the transformation matrix using the IRIS’s Geometry 
Engines. The transformed coordinates are removed from the matrix stack and the 


process is repeated until all vertices have been transformed. 


B. DOUBLY LINKED LIST 

A data structure for each vertex of the polygon is maintained in our system. 
This structure contains the vertex’s index into the vertex array. As actual 
coordinates are needed, the vertex’s index is used to extract the x, y, and z 
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coordinates from this array. A "used" flag indicates the vertexs use in a 
developed convex polygon. with a slope and y-intercept stored for the line 
segment the vertex forms with its predecessor in the array. A forward pointer 
connects the structure with the next vertex in clockwise order. while a back 
pointer connects the structure with its preceding vertex. This produces a doubly 
linked circular list of vertex structures that permits forward clockwise traversal of 


the polygon and rearward counterclockwise traversal (Figure 6.1). 


C. “WERTICES STACK 

To manage the examination and selection of vertices for a developing convex 
polygon. a stack is utilized. The size of the vertices stack is determined at run 
time when contiguous memorv is dynamically allocated. based on the number of 
vertices in the concave polygon. As a vertex 1s examined. its index is placed on 
top of the stack along with the slope and y-intercept of the line segment it creates 
with the vertex below it. If the vertex passes all selection tests. it remains on the 
stack and a new vertex is pushed on top. On the other hand, when a vertex fails 
a selection test, it is popped off the stack and the next vertex to be examined is 
placed on top. Once the selection of vertices for a convex polygon is completed, 
the vertices are passed to the IRIS system hardware and the stack 1s cleared to 


continue with the decomposition process. 
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(a) Doubly linked circular list of structures, 


Rearward us F Forward 
Pounwver Bop ee cep Pointer 





(b) Vertex structure 


Figure 6.1 
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VII. LIMITATIONS AND CONCLUSIONS 


This study has presented an algorithm for the decomposition of an arbitrary, 
planar polygon into a set of convex polygons, thus eliminating the need to 
manually perform such decomposition. With the algorithm’s dynamic allocation 
of memory. the size polygon that can be filled is limited only by the amount of 
memory available for use. The algorithm indirectly provides an extended 
capability over the IRIS system's limit of 384 vertices in one polygon (Ref. 6: p. 
zi 

However due to the time required to test each vertex for concavity, the 
algorithm is not suggested for use where real-time response is important. To 
demonstrate the overhead involved in testing for concavity. a convex polygon was 
filled and displayed 500 times by the algorithm. The same polygon was then 
filled utilizing only the IRIS’s polygon fill hardware. The results are provided in 
Figure 7.1 for various polygon sizes. 

Due to the three primary steps necessary for decomposition and the tests 
conducted on each vertex. the speed of the algorithm depends entirely upon the 
number of vertices and concave angles of the polygon. As such, the success of the 
algorithm's use in real-time display and animation is based upon a balance 
between the complexity of the concave polygons used and the speed requirements. 
We suggest that the created concave polygons be computed once and saved. 
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Figure 7.1 -— Time required to fill and 
display 500 convex polygons with 


the number of vertices indicated 
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For IRIS usage. this is perhaps best accomplished through the use o: 


mechanism of the IRIS graphics library. 
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APPENDIX — ALGORITHM SOURCE CODE 


dO doko Gk kkk ok kok ka kk kk kkk ak ag ok ak 


globals.h 


yO OK Kk kkk kk kkk kok kak kkk kk kaka kk kkk kok ak ke ok ak 


# include <g].h> 

# include <device.h> 
+ include <stdio.h> 

# include <math.h> 


[EE EE EE KOK KE KE oR eo Ek a ok a ok REE ERE ok KK KE oR EK KK KK KF 


| GLOBAL CONSTANTS 


2K 2K oo ke KK oR OK KEK KOK ok 2K OK Ek kK Kk EK OK oko Ko ok Kk ok OK KE KK ok OK KR Kk KK OK OK 
/ 


# define NULL O 
+ define BOTTOM 0 


+ define new(x) (x *) malloc(sizeof(x)) 


(6 OK Eo KE KE ok OK OK eo ROK OK ER kK KK KE KK KE oe oe ok oe ok coe EK ok ok ok KK OK kK OK 


| GLOBAL TY PEs 


2K KOK eK KK KK Ko OK oo KK RK EK ok OK RE EK KK EK KK KEK OE OK KK OK KOK OE KKK OK KOK OK A 
f 


ty pedef int boolean: 


* Define a structure to contain data about a vertex of the polygon. 
These structures will be used to create a linked list. * | 


typedef struct point { 
struct point “back; 
int index; 
float slope: 
Coord intercept: 
int direction; 
boolean failed; 
boolean segment used; 
struct point *forward; 


typedef struct point point; 
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/* Define a structure for storing data about each vertex. This will 
be used to create a stack for later use. * / 


typedef struct stack { 
int index; 
float slope; 
Coord intercept: 
int direction; 


}; 
typedef struct stack Stack: 


| 2K 2K 2 OK OK OK OK EK EE EK OK KOE OK oe OE EK EK 2K OK OK EK EE EK RE OK EK ok KEK OR OE OK OK OR OK Ko ok KE OK OE OK OK 


GLOBAL DATA STRUCTURES AND VARIABLES 


I EE EEE EE EEE EE EE EK / 


Stack *stack: 
point *start ptr. “current ptr: 
int TOP: 


boolean clockwise error flag; 
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[AIO IOI IO Eso ak doi ae deo ak ak ak ae ak 


ROUTINE: converte 
FOO COORG CO oo dodo zo aa aK aK / 
This routine takes th. initial array of vertices, passed by 
the user. and makes the calls to other routines to perform the 
necessary translation and rotations to transform the array into 
the x-y-O plane This new array is then used by the decomoser to 


breakdown the polygon into convex polygons. 


3K EK KK ERK ok EK EK KK aK aK AK RE ok ok oko oR ok OR ak oko KK aK KK KG aK OK ok 2k a KK KOK KKK KK aK / 
/ 


cpolf(total pts. pts) 


int total pts: 


Coord pts||'3); 
{ 
Int«ptl ot2, piascnt, 
float A, B. C, D, anglex. angley. anglez. *new pts: 


Matrix transforml; 
Matrix points: 


stack=(Stack *) calloc(total pts + 1. sizeof(Stack)): 


[KK EK OE EO EK EK EE KK KK OK RK EEK EK EEK KK EK OK KK 2 EK KKK OK ROKK OK KOK OK ok 


dynamically allocate memory to create an array to hold the 


vertices of the newly transformed polygon. 
KKK KK 2K KOK OK KOK KKK KK KK KK KOK OK ok KE OK KK OK OK KOK 2K KOK KK KK KK OOK KK 2K OK KOK OK OK OK OK KOK OK KOK OK KKK / 


new pts= (float *) calloc (3*total pts, sizeof(float) ); 


[REESE ER ALES EE ee ee eee 


get 3 non-colinear vertices from the polygon to be used for 


calculating the normal to the polygon. 
TK OK IK KE EK KE 3 ok KK OK KK 3K ok KK IK OK KO OE ok 3K OK KOK IK OK KK KK EK KK EK KK KK KK OK OK KOK KOK KK OK / 


extract 3 pts(total pts, pts, &ptl, &pt2, &pt3): 


[KKK OK KKK KKK KK ok ok ok oo ok KK KK KK KK KK OK KK ok ok KKK KOK KOK KOK ok KK KK OK OK OK KK OK OK ok ok ok ok Ok 


calculate the coefficients for the equations to the plane 


of the polygon and the normal to that plane. 
EAE EREERE ARERR EERE AEE EE ER EE EEE ee ee ee 
f 


determine coeffs(total pts, pts, ptl, pt2, pt3, &A, &B, &C, &D): 
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Fa Aig Re ta rt a ae ie ee oe oe Oe oe Oe ee EE 


check to ensure all vertices of the polygon lie in the same 


plane. If so then calculate the necessary angles of rotation 
ee ee rn et ee Ee A EK 


if(planar polygon(total pts, pts, A, B. C, D)) { 


FACOG IC aii II i a a 2 2k IC CIE aE a a ak I ak 


calculate the rotation about the vy axis necessary to 


bring the normal into the y-z plane. 
Saget ct ke oe eR ee Ee OE KE oe 


angles from normal(A, B, C, ’y’. &anglex. &angley); 


nei ee ee ee ee EEE EEE 


create the transformation matrix necessary to multiply 


with the array of vertices to generate the rotation. 
Saar RE EE EE OE a 


build transform matrix(transform], angley, ‘y’, pts, pt1l); 


Ri ee ee Ee Ee Ee ee Ee 


perform the matrix multiplication with the array of 
vertices and store the new coordinates in the new array 


new pts. 
oa Et Ee EE eR EE 


map pts(total pts, pts, transform1, points, new pts); 


ee SE ES EE EE EE EE 


now extract three non-colinear vertices from the new 


array of vertices, new pts. 
FEES EI a a a aaa tallied | 


extract 3 pts(total pts, new pts, &ptl, &pt2, &pt3); 


PEE EAA EE EE EEE EEE A OEE EE HE ee OG EE ke Se ee ok ok a oe 


since the polygon has been rotated about the y axis we 
calculate the coefficients for the equations to the 


polygon one more time. 
NN IIIS SEDATE ta aladdin iadiadadalaldataladall | 


determine coeffs(total pts, new pts, ptl, pt2, pt3, &A, &B, &C, &D); 
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ok eK Ko KK KK KE KE EK KK EE OK KE ee eR ee ie 


calculate the rotation about the x axis 


Dd ia dK Kk dR kok ack ak daa ok ae ka 
i 


angles from normal(A. B. C. °*x’. &anglex, &angley ): 


2K KKK ok KK OK ok RK Kook oo KOK KK ok Kok ok OK OK ok ok 2 KOK ook OK OK ok ok ok OOK ok oko KK ok ok OK OK KOK OK ok & 


create the transformation matrix 
Te OK SOK OK OK KOK KK OK OK KK KE ok EE OK EK ok ROK OK OR OR KE ok ok eK OR SK OR oe es kk ee ak er 


build transform matrix(transform1. anglex, ’x’. new pts, ptl); 


[8 EE Eo oo eo EE EE Oo ERE aE EEE ook oo KE ORK CK KO ok ok ok kK 


perform the matrix multiplication 
RO RK RR EK OK OK KO oe KOK KK ok KOK OR ok OKO KOK KOK OK KOK KK OR ok KOK KOK OK OO KOK OK ok OK OOK OK ok ok KK KOK KK 


map pts(total pts, new pts. transform]. points. new pts): 


RK OK OK RK ROK OK KK ok KE kK KK KK EK KKK OK KK OK OK ok OK OK OOK OK KK KK KOK KK KE 


call the routine which will perform the decomposition 
upon the newly transformed polygon which now lies in the 
x-y-0 plane.(thus to be handled as a 2D polygon) 


OK KE KK ok KKK 2 oe KK OK ok OK ok OK Kok Ko OK OK KOK KK ok OK OR EO KK KOK KOK KOK ok OK OR Kk OK OK KKK KK KKK 


split(total pts, new pts, pts): 


Ae KK Eo OK KK EK ook oo KK KK aK ok EK RK ok oO KK OK ok OK ok KOK ok OR 2K KK 2K OK OK OK EK ok OK OK ok KOK KOK 


once the dynamically allocated memory for the array 


new-pts has served its purpose. the memory is freed. 
ho KKK ok OK oe ok Kok 2 oko Eo ok KO KK 2K OK ok ok KOE ok ok OK ok Rook oo Ko oo ok oo Ko Kok KKK KK KK / 


cfree(new pts. 3*total pts, sizeof(float)); 
cfree(stack, total pts, sizeof(Stack)); 


3K OK ok Ko KK KK ok OK OE ok ok OK OE ok oo eK Ko EK ok ok ok KE ok oe 3 KK oO KK KK OK ok OK OE KOR OK OK KOK OK KOK OK Kk 


if not planar then the decomposition is not attempted. 
oh KOK ok ok ok ok oko 3K OK ok ok Ko oe oo kook EK Ko oO KK 2K ok ok ok KK KK EK ok ok OK KOK oe OK oo ok OK ok Kok KOK OK OK OK 


return; 
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[2 EE EK EE EE EE EK IE 2 EER EOE EEE EO EE EE EE ER ao EE EK KOK KF 


BROW TING: extractc 


OK KOE KK OK OK EE EEE KE EEE EE RE EE OE EE EE ERE KK KI KE EE KK 


This routine searches the array of vertices to find three points which 


do not lie on the same line and returns the indices of these points 
a a EK EEE EE EEE EEE EERE 


extract 3 pts(total pts, pts, ptl, pt2. pt3) 
int total pts; 


Coord pts|][3); 
ime ptt, pt2, “pts; 


{ 


int count, finished: 
Soei— 0) pi2— 0. pts— 0; count— 0: finished— 0: 
/* The first vertex in the array (index=0) is selected for the first point. */ 


while(count < total pts && !finished) { 


/* We then search the rest of the array of vertices to locate two more * / 


count+= lI: 
if(*pt2 == 0) { 


/* If the second point has not been selected yet then we check that 
x. y. and z coordinates, of the vertex currently under exam, are 


different. * / 


if(pts|*pt1]/0! != ptsicount]{0] || 
pts *ptl]/1) != pts/count}]|1] || 
pts'*pt1][2] != pts/count][2] ) 
/* If so then the vertex becomes the second point. */ 
*pt2= count; } 


else if(*pt3 == 0) { 


/* If the third point has not yet been selected, then we check that 
the coordinates of the vertex under examination are different 
from both the first and second point. */ 


if((pts|*pt1]/0) != pts/count] 0} || 
pts|*pt1)/1] != pts[count][1] || 


39 


pts|*pt1]/2) != ptsicount|/2. ) && 

(pts *pt2] 0 != ptsicount! 0, 

pts(*pt2][1] != ptsicount]|1) | 

pts|*pt2|/2 != ptsicount||2) )) 

/* If the coordinates are different then the vertex is selected 
as the third point. */ 

*pt3= count: } 


else 


finished= 1: 


return: 
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[KE EEE KK EK EE EE KK EEK EE EE EE EE FE EE EE EK 


ROUTINE: is planar.c 


OO ao OK dk cde kkk kek / 


This routine takes each vertex of the polygon and uses it to solve the 
equation for the plane as calculated earlier. If every vertex lies 


in the plane. then a value of 1 is returned. 
2K eK 2 EK 2 EE EE EE EE ee EK ee 3 Ee SC ee EE EE EE ee 2 EC ee oe oe EE EE / 


planar polygon(total pts, pts. A. B. C, D) 
int total pts: 


Coord pts|j/3); 
float A, B.C, D; 


int count, planar: 
count= Q; planar= 1: 
while(count < total pts && planar) { 
‘* The equation for the plane is setup to equal zero if the coordinates 
of the vertex lies in the plane. A range of .02 on either size of 
zero 1s allowed to handle round off error. This test is done for 


each vertex in the polygon. If any vertex fails the test then the 
function returns (0) for non-planar. */ 


if((A*(pts count|/0/)+ B*(pts/count|/1 )+ C*(pts:count]|2))- D) > .02 || 
(A*(ptsjcount| 0})+ B*(pts)count]/1!)+ C*(pts count]{2))- D) < -.02) 


planar= 0; 


count+= I: 


} 


return( planar); 
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[KOK IK OK KOK OK KK OK KK KK 2K 2K OK KK OK KK 2K KEE OK KK oR OK ROK OK OR KK KOK OK OK OK OK KOR KOR OK KOK KKK KKK 


ROUTINE: coeffs.c 

BO OOO oo kk ok OG ok Oo kk a ke 
This routine uses the three vertices selected from the polygon to 
calculate the coefficients that will be used in both the equation 


of the plane the polygon lies in and the normal to that plane. 
SE SE EE Eo aE EE EE EK EO EE EO I EE EE EE EE EO EE EE EE EE EO Ek KEK / 


determine coeffs(total pts, pts, ptl, pt2, pts, A, B, C, D) 


int total pts: 
Coord pts/||3}; 

int ptl, pt2, pt3; 
fleat Aja B. Ce Db: 


{ 


Hoat x0. x1y x2. 0e yey 270. 72 


x0= pts ptl |0.: 


yO= pts pt1]/1/; 
z0= pts(/pt1]/2]; 


x1= pts pt2!/0:; 

y l= pts|pt2]|1; 

z1= pts|pt2)|2); 

x2= pts /pt3]|0;; 

y2= pts|pt3]|1: 

z2= pts|pt3]/2|; 

*A= (yO-y1)*(z2-z1) - (y2-y1)*(z0-z1); 

*B= (20-21)*(x2-x1) - (z2-z1)*(x0-x1); 

*C= (x0-x1)*(y2-y1) - (x2-x1)*(y0-y1); 

*D= (*A)*pts[pt1](0} + (*B)*pts/pt1][1] + (*C)*pts|pt1)[2); 


Fevury, 
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DK OREO OKOKER SOS kkk ok ok ok 


ROUTINE: angles.c 

2K KK KK KK KK KK KK ok KK ok KK KK KK ok OK ok OK KK KE Ko ok KK oR KOK OK OK OK OK KKK KKK KOK KK ok ok ok ok ok OK 
This routine utilizes the coefficients. determined from the normal to 
the polygon. to determine the angles of rotation necessary about the x 
and y axes. These rotations are necessary to bring tne polygon into the 
x-y plane. This routine is called the first time to calculate the 
rotation about the y axis. The second time it 1s called is to determine 
the x axis rotation. The parameter "axis" is used to pass which angle 


is to be calculated. 


oe EE KOE ok EK KE kK OE KKK IE EEE EK KE EK KK EK KK EE KKK EK ok oe ok Ke OK OE OK oe 
/ 


angles from normal(A, B. C, axis, anglex, angley) 


float A. B. C: 
char axis: 
float *anglex. *angley: 


{ 


float pi, degx. degy. degz: 
pi= 3.1415926536; 


degx= acos(A/sqrt(A*A+B*B+C*C))*360/(2* pi): 
degy= acos(B/sqrt(A*A+B*B+C*C))*360/(2* pi); 
degz= acos(C/sqrt(A*A+B*B+C*C))*360/(2* pi): 


if(A==0 && B==0) { 
/* simply translate the polygon to the origin * 
*anglex= 0.0: 
*angley= 0.0; } 


/ 


else if(A==0 && C==0) { 
/* rotate polygon 90 degrees about x axis */ 
*anglex= 90.0; 
*angley= 0.0; } 


else if(B==0 && C==0) { 
/* rotate polygon 90 degrees about y axis */ 
*anglex= 0.0; 
*angley= 90.0; } 


else { 
/* must calculate amount of rotation about the x & y axes */ 
“anglex= acos(C/sqrt(B*B+C*C))*360/(2*pi): 
*angley= acos(C/sqrt(A*A—C*C))*360/(2*pi); 
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if(axis == ’y’) { 

if(degx < 90.0) 
*angley= -(*angley): 

*anglex= 0.0: } 

else { 
if(degy > 90.0) { 
“anglex= -(*anglex): 
*angley= 0.0; } 

} 

} 


/* convert the angles into tenths of degrees as needed by the IRIS system. 


*anglex= 10*(*anglex); 
*angley= 10*(*angley): 


return: 


44 


x 
} 


ee ee Pe ee ee a a ee ee ee eS oe 


ROUTINE: transmatx.c 


Se ee ae a ge shee seo deed ston ance ae oe Se Ged oe oe ec ae de ae eee Se ae eon on ae eae oa 


This routine takes a matrix structure and initializes it to an identity 
matrix. Then the necessary translation values and rotation Values are 


added to the matrix to create the resultant transformation matrix. 
FCC CIC CI ad a I I I ICICI Ea / 


build transform matrix(transforml, angle, axis, pts, pt1l) 


Matrix transform1: 
float angle: 

char axis; 

Coord pts|][3,; 
mt) ptt’: 


\ 
it 1, ): 
/* create the identity matrix *| 
for(i=O; 1<4: i~=1) { 
for(j=0: j<4; j+=1) { 
iiti=—)) 
transform] i] j = 1: 
else 
transform] i||} = 0; 
} 


a i 
j 
pushmatrix(): 
loadmatrix(transform 1); 


/* add in the rotation values */ 
rotate((int)angle. axis); 


/* add in the translation values */ 
translate(-pts|pt1]|/0), -pts|pt1]/1), -pts/pt1) 2]); 


getmatrix(transform1); 
popmatrix(); 


return: 


} 
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JIC ICICI IORI ook ck tok ek sek ae ak 


ROUTINE: mapover.c 

FOO OOO III IOI ORI HIGGS ICI rac or ack acca eae teak ak / 
This routine places up to 4 vertices into a matrix at one time then 
performs a matrix multiplication with the transformation matrix to 
produce x, y. and z coordinates for the polygon that is translated 
and rotated some amount about the x and y axes. These new coordinates 


are restored in the array new pts for further manipulation. 
hE GOO IOI GOK Oa kok ko ak ag kag ak ak ok 


map pt-‘total pts, pts, transforml, points, new pts) 


int total pts: 

Coord pts|||3; 

Matrix transform|, points: 
float new pts|! 3 ; 


int count. count2, cnt: 

count= 0: count2= 0: 

while(count < total pts) { 
cnt= OQ; 


while({cnt < 4) { 
/* place the coordinates of a vertex into the matrix */ 
points|cnt||O| = pts count] 0): 
points|cnt|/1) = pts count] 1): 
pointsicnt]|2) = pts:count}|2); 
pointsicnt|/3) = 1.0: 


if({count+= 1) < total pts) 
cnt-—= |; 
else 


break; 


/* save the state of the matrix stack and perform the multiplication */ 
pushmatrix(); 


loadmatrix(transform1); 
multmatrix (points); 


getmatrix(points): 
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popmatrix(); 
cnt= 0: 


while(cnt < 4) { 
/* move the newly calculated coordinates back into the vertices array */ 


new ptsicount2]'0/= pointsicnt|/0: 


new ptsicount2| 1)= pointsicnt|/1); 
new pts'count2]\2'= pointsicnt|/2 ; 





if((count2—-= 1) < count) 
cnt+= 1; 
else 


break: 


return; 
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15 oR RRR SE EO HCO SER oi 9 2 Ee ke eR a ae ee cee 


ROUTINE: split.c 
SOKO Soo ok ok ok kkk kkk kak: / 


é 


This routine takes the array for the transformed polygon and calls 
the necessary fu tions and routines to decompose the polygon 
and pass the vertices of the convex polygons to the IRIS system 


for filling. 
KER KKK KKK KK KK KKKKKKEKEKKAKEKKKKKKKKKKKKKKKK EK Ke a a 


split(total pts. pts. old pts) 


int total pts: 
Coord pts!/]|3); 


Coord old pts/|'3); 


Object decompose polygon(): 
Object whole thing: 

boolean clockwise: 

point *pl, *p2: 


SERRE E LAER RE EAR EEK EA EEE EEK EKER KEE EEE EEE EEE 


Create the doubly linked circular list. "main chain", for the 
concave polygon represented by the pts array. Also initialize 


flags and ptrs that will be used in decomposing the polygon. 
OK OK KKK OK KOK Ko KK KK KK KK OK ok ok KKK KK KK oR KOK KKK ok OK ROK KKK KOK OK OK ROE KK EK KK KKK / 


é 


clockwise= clockwise ordering(total pts, pts); 


create main chain(clockwise, total pts, pts); 


[KK KE KK OR EK KE KK KK KK a EE oO KK KE OK OK KOK OK KE OK OK ORK KOK KOK OK KOK KO 


Decompose the main polygon and build the set of simple convex 


polygons. 
TE EE KE RR OK EEK KO EE OK KEE KK ok oe KK KO KK KR KK RK RR OK aR KOR KK RK OK KO / 


‘ 


decompose polygon(total pts, pts, old pts); 


return; 
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DOC Colcord cori dco rE dered de ara dak a ae a 


ROUTINE: main chain.c 


ok OOO oko oi ooo a dock ic kak ak ie ak oe / 


f 


Create the doubly linked circular list. "main chain". for the 
concave polygon represented by the pts array. Also initialize 


flags and ptrs that will be used in decomposing the polygon. 
Te KK KK Ok OK 2K ok Kok Ko KK OK ok ok eK oko oo KK oo oo oi oi oo EE OK CE OE EE ok Eo ok oi IE oe EE Eo / 


create main chain(clockwise. total pts. pts) 


boolean clockwise: 
int total pts: 
Coord pts|//3;; 


int countl: 
Hone pil, p2: 
/* allocate memory for a structure to contain data on the first 
vertes of the polygon * 


pl= new(point); 
/* set pointers to this structure and initialize all values */ 


Start ptr= pl; 

current ptr= pl; 
pl->index= 0; 
pl->direction= 0; 
pl->failed= FALSE; 
pl->segment used= FALSE: 


/* if polygon vertices are ordered clockwise then we index the vertices 
array beginning at 0 and incrementing from there. Other we must begin 
at the other end of the vertices array to extract our coordinates. In 
this way a forward move along our doubly linked list will equal a 
clockwise move around the polygon. */ 


if(clockwise) 
countl= |; 
else 


countl= total pts- 1; 
/* for each vertex of the polygon we build a point structure, calculate 


the slope, intercept. initilize fags and connect the structure 
via pointers to its preceeding and succeeding vertices in the 
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polygon. The result is a doubly linked list of structures. * 


do{ 
p2= new(point); 
pl->forward= p2; 
p2->back= pl: 
pl= p2: 
pl->index= countl: 
determine slope(pts pl->back->index)'0), pts|p1->back->index||{1,, 
pts pl->index)}0), pts pl->index]/1), &(pl->slope), 
&(pl->intercept). &(pl->direction)): 
pl->failed= FALSE: 
pl->segment used= FALSE; 
if(clockwise) 
countl+= 1: 
else 
countl-= 1: 
}while(count]< total pts && countl > 0); 
‘* connect the first vertex structure to the last vertex of the polygon 
and calculate the slope and intercept of the line segment the two 


create. * 


pl->forward= start ptr; 

Start ptr->back= pl; 

pl= start ptr; 

determine slope(pts pl->back->index]'0,, pts|pl->back->index)/1 , 
- pts'p1->index) 0). pts pl->index|/1), &(p1->slope). 
&(pl->intercept), &(pl->direction)); 


return: 
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[KKK KK OK KK KK OK OK KOK OK ok OK OK ok ok OK KOK OK OE KK OK OK OK EK Ko OK ok OK OK ok OK OE OK OK OK OK OK OK ok OK ok OOK OK ok OK OK KOK oF 


ROUTINE: clockwis.c 


EK KK KK KK KK ok KK KKK ok ok ok ok ok IE OK KE aK OK aK OK OK ok ok ok KK oe Ko 2K aK OE OK OK OK ok OE KK i OK OK ok KK KK / 


This routine calculates the area of the arbitrary planar 
polygon once it has been transformed to the 2D x-y plane. 

A negative area indicates a clockwise ordering of the vertices 
while a positive area indicatesa a negative orderig of vertices 
A "0" is returned for counterclockwise ordering while a "1" is 


returned for clockwise ordering. 
2K KK 2K IK ok IE KK EK KK KE OK kK ok a KE oe KE KK ok ok ok KK KE KK eK KR aK KK OK kK aR KK ok ok / 


clockwise ordering(total pts, pts) 


int total pts; 
Coord pts|||3); 


float area: 
int 1. result: 


area— 0.0; 


for(i=O: i< total pts- 1; i+= 1) { 
area= area + pts|i] 0 * pts|i+1][(1); 


} 
area= area + pts total pts- 1/0 * pts/0//1: 


for(i=0; i< total pts- 1;i+= 1) { 


area= area - pts|i] 1 * ptsji+1,|0]; 


} 
area= area - pts|total pts- 1]/1] * pts[0|'0); 


if(area >= 0) 


result= 0; 
else 
result= 1; 


return (result); 


a1 


[KOK KKK ROOK ROKK KOK KK KORO OK Kok ak KOK KK kK Kk aK kK kkk kkk 


ROUTINE: decompose.c 

OOK kok kkk oO KK kk kok ko KK a kkk / 
Decompose the main polygon and as convex polygons are determined 
these are passed through calls to the IRIS system to be filled. 


SOK KKK KKK OK KKK Kokko KKK KK Kk RK OK KKK ak kK kok bok ok ok ok kkk ok kok kok / 


Object decompose polygon(total pts, pts. old pts) 


int total pts: 
Coord pts|].3); 
Coord old _ pts|/[3': 


boolean no point found. search incomplete. adjust needed. Ok; 
boolean skip flag. skip search: 

int pti. pt2. poly count. value: 

point “working ptr: 


clear stack(): 
get pt(total pts, pts): 
get pt(total pts, pts): 


search incomplete= TRUE; 
skip search= TRUE: 
poly count= 0: 


fo EK KK EK OK OK OK OK 2 OK OK 2K OK OK 2K 2K 2K OK OK OK KOK KK OK KK OK OK OK OK KK OK KK OK KK OK OK KOK KK OK KK OK OK OK K 


Once an initial two points are taken from the concave polygon. 
the process moves through the main chain to select out points 
for constructing the simple convex polygons. Calls are made to 
various routines to check for validity of pts that are accepted 


into a polygon. 
ne KK KOK KK KOK KK OK OK OK ok KOK OK OK OK 2k OK ok OK 3K KOK OK OOK OK OK KOK OK KOK OK OK OR KOK OK OK OOK KOK OK OOK OK OK KOK KK KK KK KKK / 


while(search incomplete) { 
adjust needed= TRUE: 
OK= TRUE; 
while(OK) { 
get pt(total pts, pts); 


/* we first check to see if we have made a complete loop of the 


vertices. */ 


if(loop completed()) { 
working ptr= start ptr: 
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/* find the structure in the linked list that matches the 
vertex on top of the stack *. 


while(working ptr->index != stack TOP!.index) { 
working ptr= working ptr->forward; 


j 


/* check to see if vertex creates a line segment with the 
vertex below it. on the stack, that intersects any edges 
of the concave polygon *, 


if(stack TOP-1).index != working ptr->back->index) { 
value= intersecting lines(ptsistack’ TOP|.index)'0), 
ptsistack TOP .index| 1), pts/stack/TOP-1 .index|/0,, 
pts|stack! TOP-1).index'!1), stack; TOP! .slope. 
stack TOP .intercept. total pts, pts): 


if(value) { 
/* if vertex creates an intersecting line segment then we 


remove it from the stack and reset the pointer to the 
linked list * ’ 


current ptr= start ptr; 
while(current ptr->index != stack/TOP-1).index) { 
current ptr= current ptr->forward: 


} 


current ptr= current ptr- >forward; 
release pt(): 
release pt(); } 


else { 


/* since the line segment create is OK, we simply remove the 
vertex from the stack since it is the same as the one on 
the bottom of the stack, namely the beginning vertex */ 


OK= FALSE; 
release pt(); 


} 
} 


/* if we have not completed the loop then we must conduct the 
tests to ensure valid selection of vertices for a convex 


a3 


polygon * 


else if(pt out of bounds(totai pts, pts)) { 
release pt(); 


} 
} 


/* if at least three points have been selected then... * 


if(valid pts()) { 
/* save the beginning and ending vertices for future reference * / 


ptl= stack BOTTOM index; 
pt2= stack TOP .index: 


‘“ pass the vertices to IRIS system hardware for filling and display 


define polygon(old_ pts): 
poly count—= 1: 


-* if this is a normal decomposition and not a "fill in" polygon to 
cover holes left in the concave polygon the the do the following *. 


if(skip search) { 
‘/* set the flags that indicate next time through this loop we 
will be working on a "fill in" polygon */ 


adjust needed= FALSE; 
skip search= FALSE: 

‘* adjust the pointers to the linked list so next time through 
the loop we have the start and end vertices of the previous 
polygon to begin working with */ 


working ptr= start ptr; 
while(working ptr->index != pt2) { 
working ptr= working ptr->forward; 


} 


/* if the start and end vertices of the previous polygon are 
adjacent on the concave polygon then forget about building 
a "fill in" polygon */ 


if(working ptr->forward->index != pt1) { 
while(current ptr->index != ptl) { 
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Bk ] 


current ptr= current ptr->forward: 
} 
get pt(total pts. pts): 
current. ptr= working ptr: 

get pt(total pts, pts): } 

else { 

adjust needed= TRUE: 

} 

i 
} 


/* if the polygon is not fully decomposed yet we select 
two new vertices to continue the decomposition process * / 


if(adjust needed) { 

skip search= TRUE: 

if(object complete(total pts)) 
search incomplete= FALSE: 

else { 
adjust search area(pts): 
clear stack(): 
get pt(total pts, pts): 
get pt(total pts. pts): 


revurn: 
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DOK OG KICKER GR kK dk kk ak 


ROUTINE: clear stack.c 


kK EK Gk Kk kkk kkk dca ke / 
/ 


This routine simply resets the stack pointer to -1. In this 
manner when the stack has its first item pushed on, the pointer 


will accurately reflect a value of 0. (BOTTOM OF STACK) 


KKK KR KKK 2K Kk KK KK KK OE 2 2 2K Kk ROK ok OE EO KOK KKK KOK KKK KK KK KEKE KK KK KK KK KKK / 
f 


i 


clear stack() 


TOP= -1: 
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[ROK KOK KK ok ok ok OK OK EE EK KE EO ROK OK OK OR KE oR OK OE EO OK OOK OF ER OK KK ROK Ok ok 


ROUTINE: get pt.c 


ee eee ee eee eee KE ee ee a Ke / 


This retrieves the next point in the main chain. as identified 

by the current ptr. and places the point on the top of the stack 
In doing so it increments the TOP variable. Also as the point 
is retrieved. the slope. intercept. and direction of the line 
segment it creates with the vertex below it on the stack is 


calculated. 
ko OK OR kok ok kok oK a kokok kK dokokak dokok ok kk ok kkk / 


get pt(total pts. pts) 


int total pts: 
Coord pts/|/3 ; 


{ 
TOP-+= 1: 


if(TOP <= total pts) { 


'* select the vertex in the linked list that is currently pointed to *, 
stack TOP .index= current ptr->index: 

/* and then advance the list pointer *, 

current ptr= current ptr->forward; 

nm DOP > 0) { 

/* take the x. y. and z coordinates for both the vertex on top of 
the stack and the vertex just below it. to calculate the slope 
and y-intercept of the line segment they create */ 

determine slope(ptsistack TOP-1).index||0). pts\stack/TOP-1’.index||1,, 

ptsistack TOP. .index||0,, pts'stack! TOP].index|/1). 
&stack/ TOP|.slope. &stack TOP |.intercept, 
&stack/TOP).direction); 


} 
} 


else { 
printf("STACK overflow occured while retrieving another point"); 


} 


return: 
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DO OK OKO KOK OK KOKO KKK RK KKK aK kkk kk kk de ok 


ROUTINE: release pt.c 


OOK KCK OK OK OK OK KK OK KK OK OR IKK Kk kK kK doko / 
/ 


This routine removes the item on top of the stack by 


decrementing the top of stack pointer TOP. 


26K KK OK ORK KOK OK KOK KOK OK OK Oe OK OK Ke KK OK ok KE OK OK KKK KOK 2K OK OK KOK OK ROOK OK KK OK OK KOK OK OK KKK KE KKK KK / 
/ 


release pt() 


{ 


if(TOP >= BOTTOM) 
TOP-= 1: 

else { 
printf("STACK underflow occured while trying to release a non-"); 
printf("existence point."): 


} 


return: 
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1K OK KO OK KOK OE KE OK KK ok ok ok KK ok OK OK KOK KOK KKK KOK KK KK KOK SK OK OOK OK OK KK KR KOK KOK OK OK KOK OK OK OK KK 


ROUTINE: slopes.c 
EP RELA Ee ER EE He EY 
This works through the pts array. and the stack array. By 
taking value of the TOP of stack and the point beneath it. 
the corresponding coordinates are obtained from the pts array 
and a slope, intercept and direction are calculated. These 


values are then stored in the stack. 
a lara ERE ER EER EEK E 


determine slope(xl, yl. x2, y2. slope. intercept. direction) 


@oord xl, vl. x2, y2: 
float “slope: 

Coord “intercept: 
int *direction: 


Coord dx, dy: 


/* calculate the delta x and delta y for the to points * 


dx= x2- x1; 
dy= y2- yl; 


/* if the line segment is vertical then use the value 99999.0 to 
indicate such (avoiding division by zero problems) * 


if(dx == 0) { 
*slope= 99999.0: 
*intercept= x2; } 
else { 
*slope= dy/dx: 
*intercept= y2 - (dy/dx) * x2: 


} 


/* based upon the changes in delta x and delta y we determine if we are 
traversing down a slope or up a slope */ 


Miay < 0) { 
if(dx < 0) 
*direction= 1; 
else 
*direction= 0; 
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else { 
itax 290) 
*direction= 0; 
else 
*direction= 1: 
} 


return; 





koko ok OK kkk kkk kok ok okok kok kkk ko kk kok kok kok bk ko kk dk kk ak kde kot ae 


ROUTINE: compare.c 
Be ee oe ee ee ee Ee ok SE EO OE OE EE Ee OO OO 2 OE SE OK oo EK 2 Ko 
This routine determines if a point lies to the left or right of 
a line segment. This is done by using the x and y coordinates 
of the point with the slope. intercept. and direction of the 
line segment to create a parallel line segment. The y intercept 
(or in the case of vertical line, the x intercept) is used in 
conjunction with the direction of traversal along the line 
segment to find which side of the line segment, the point lies 


on. 
ok KOK EK KKK KK Ko oo SE Ko ok ao RK ok ok ok ok ok aR ok ok oR Ko ok oo oR ER aK a EEE ok oo 


compare boundry(x1.y1.slope.intercept.direction ) 


©oord xl. yl: 
float slope: 
Coord intercept: 
int’ direction: 


float test value: 
int value: 


/* if the line segment is vertical then we use the x-intercept to 
determine if a vertex lies to the left or right of another line 
segment */ 


if(slope == 99999.0) 
test value= x1; 
else 


test value= yl- slope * 


Sols 
/* based upon traversal up or down a sloping line we determine if 
a vertex is to the left or right of the line segment */ 


switch (direction) { 


case 0: if(test_ value > intercept) 
value= I: 
else 
value= 0; 


break: 


case 1: if(test_ value < intercept) 
value= 1; 
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else 
value= 0; 
break: 

} 


return (value): 





OK OK OK ROR OK kok ok aR kk cok ae sk aoe oko ak ook de ok ak de ok ak ke ak 


ROUTINE: pt OOB.c 


SOKO ROR ORK ROR RRR RS RK ER RE Rk ok kk kak a a ak 2 


This routine takes the vertex index from the top of the stack 
and uses its x. y coordinates to determine if the point lies 
out of bounds with respect to the convex polygon being 


created. 
OK KOK RRO KKK RR RRR kkk oR RR kk RE ER ae a a oe oe ok ok of 


pt out of bounds(total pts, pts) 


int total pts: 
Coord pts;|/3); 


{ 


Woord xl. vl: 

int ptr. value; 
boolean not finished: 
point “working ptr: 


* obtain the coordinates for the vertex on top of the stack * 


x1 = ptsistack!TOP).index!!0'; 
vl = ptsistack/ TOP).index) 1; 
ptr = TOP- 1: 


value= 0: 


‘* use the above coordinates to determine if the vertex lies to 
the left (out of bounds) with respect to any line segment 
created thus far in the selection of vertices */ 


Misia (per = BOTTOM && Wvalue) { 
value= compare boundry(x1,y1.stack|ptr|.slope.stack|ptr|.intercept. 
stack ptr|.direction); 
ptr-= 1: 


} 


[EK EK Ee EE OE ER OO oe EK EO Ko ok OK oe Eo KK CE oo EK oe ok a I OE 2K OK ok ob 


If point is not out of bounds with any line segments created 
then we check to see if the segment produced by this new point 
addition causes the start point for the polygon to fall out of 


bounds. 
ee eg OER ER EEX f 


if(tvalue) { 
xl= pts|stack/BOTTOM|.index)/0;; 
y1= pts|stack; BOTTOM) .index|/1); 
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value= compare boundry(x1.yl.stack TOP .slope.stack| TOP! intercept. 
stack TOP direction): 

if(value) { 
release pt(): 
current ptr= current ptr->back:; 
working ptr= start ptr; 

} 

j 


3K Oh ok ok OK ob oo OK 3K 2K KK KK OK 2K KK 2K 2K OK OK 2K ok OK KK ORO KOK oR OK ok OR Ko OK oO OK KS OK OK OK ok OK ok OK oe OK KKK ok KK 


If the newly created line segment does not cause the initial pt 
of the developing polygon to fall out of bounds. then we check 
to see if the line segment intersects any line segments of the 


concave polygon. 
5h OK RR KOK OK ok OK oR a aK OK EE OK OK ok ak OK ok oe ok Ko aK OK aK a ok oo ok ok oe oe OK OK oR KKK oo ok oR OK KK OK KOK OK KOK KKK KOK OK 


if('value) { 
working ptr= start ptr: 
while(working ptr->index != stack|TOP).index) { 
working ptr= working ptr->forward; 


if(stack TOP-1'.index != working ptr->back->index) { 
‘* we use the coordinates for the vertex at the top of stack 

and the vertex below it to create a line segment which we 

then check for intersection with anv of the edges of the 


original concave polygon * 


value= intersecting lines(ptsstack' TOP .index],0 . 
pts'stack TOP).index||1 . ptsstack|TOP-1 .index!|0. 
pts stack TOP-1).index)|1), stack TOP .slope. 
stack TOP].intercept, total pts, pts); 


[OR Ro OK KK ok oR OK ok oo 2K oR ok eo OK oR ok OK ok KK KK KK OR oR KR oo oo ok ok Kk ok Ko KK OK OK OK OK KOK OK KOK 


If the newly created line segment does not intersect any line 
segments of the concave polygon then we check to see if the 
line segment of the original polygon, that the top of stack 


point heads, causes the point below top of stack to fall out. 
Be 2 2 Oe GS 2 ee a aE a a a ee Ae Ae a ee A Ae A A ee 


if(!value) { 
if(current ptr->back->back->index != stack|TOP-1].index) { 
x1= ptsistack! TOP-1].index] 0): 
y1= ptsistack! TOP-1].index]/ 1}; 
value= compare boundry(x1l,y1l.current ptr->back->slope. 
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current ptr->back->intercept. 
current ptr->back-> direction); 


} 
} 


a eee eee ee i hk eo hk eh ob a ed ok EK OK + ok 


If the point 1s determined to be out of bounds. then the 
appropriate pointers are adjusted to continue the decomposition 


process. Otherwise the pt is a valid addition to the stack 
ko KOK kk kK koko ho kokok honk a KE kok kok bok Kok kOK Rok ag ak oak nc a ok 


if(value) { 
working ptr= start ptr: 
not finished= TRUE: 
do { 
if(working ptr->index == stack!TOP .index) { 
not finished= FALSE: 
if(working ptr->failed == FALSE) 
working ptr->failed= TRUE;} 
else { 
working ptr= working ptr->forward: 
} 
} while(not finished): 
} 


return(value): 
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[KKK KK OK OK KK KK OK OK KK KOK EO KOK RK ok KK EK ROK KK KKK EE KES 


ROUTINE: loop complet.c 

3K KK 2K OK OR 3 OK KOK ok ok OK 2k OK 2k ke KK OK OK ok OK OK kK dc 2K 3k OK OK OK OK OK OR OK OOK OK kK OK KE OK OK OK EK om OK EK OE KK 
This routine checks to see if the vertices of the concave 
polygon have been completely traversed while developing a single 
convex polygon. This is done when the vertex on top of the 
stack is the same vertex on the bottom of the stack. When this 


occurs a Value of | is returned. 
SK OK SK oe SRK OK oR CE EK Ke SK SE She Se hc She SR Shh Se She heh ae hc ie ic ic ct crc iel ace ae 


int loop completed() 


{ 


int value: 


(TOPS =1) 4 
if(stack TOP .index == stack BOTTOM..index) { 
value= 1; } 
else { 
value= 0: 
j 
} 


else { 
printf("Error in stack index while checking TOP & BOTTOM values"); 
value= -1l: 


} 


return (value): 
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ROUTINE: adjust.c 
OOK ROK Ok OR kkk kok kkk kK ok ok 2 kk dk ak ok 
This routine searches the circular list of vertices. checking 
the various condition flags, to locate a vertex or line segment 
of the originial polygon which has not been utilized. As one is 
located, the list’s currency pointer is moved to that structure 


in the list so that further polygon decomposition may continue 
FOO OCR GO AOI CS CSCC CC II dR I I I a a | 


adjust search area(pts) 


Coord pts|]'3 : 


boolean done: 
pl: 


> 


point 


done= 0: 
pl= start ptr: 


do { 
if(pl-> failed == 1) 
/* check to see if the line segment to either side of the vertex 
has been used in a previously developed polygon. * 
if(!(pl-> forward-> segment used) !(pl-> segment used)) { 
current ptr= pl: 
pl-> failed= 0; 
done= 1; } 
else { 
pl= pl-> forward; 
} 
else { 
pl= pl-> forward: 
} 
} while(!done && pl != start ptr); 


return: 
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ROUTINE: valid pts.c 


en ee ee ee ee ee ee ee ee ee ee ee eee ee ee ee ee ee eee ee ee ee ee 


This routine checks to ensure that at least 3 vertices have 
been determined for the convex polygon. If not then the 


value of O is returned. 
KK OK KARE KKK KKK KOK KKK KKK KKK KK KEK KEK KKK KKK KK KKK EK KEK KKKA KKK KE KEK KEK KE / 


boolean valid pts() 


{ 


int value: 


if(TOP >= 2) 


value= IL: 
else 
value— 0: 


return (value): 
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HOO ho a dh ob ka rk de ka ae ok dee kok 


ROUTINE: define poly.c 


COORG OO Ok oo doo bod kkk kaa 4 


This routine takes the vertices accumulated on the stack and 
uses the IRIS features pmv and pdr to actually perform the 
display of the polygon created. The index for each vertex 

is used to access the original array of points to get the x. vy 


and z coordinates used. 
SOK OOK cho kk oi Khoa oka kok ack kok ok kk ak 


define polygon(pts) 
Coord pts 3: 


{ 


Int count: 
point ‘track ptr: 


track ptr= start ptr; 


/* pass the starting vertex for the newly created convex 
polygon to the IRIS hardware with a point move command * 


pmv(ptsistacki BOTTOM ..index| 0. 
pts stack BOTTOM ..index] 1). 
pts stack BOTTOM index! '2]); 


/* locate this vertex in the linked list *. 


while(track ptr->index != stack/ BOTTOM. .index) { 
track ptr= track ptr->forward: 


} 


/* set the flag to indicate the vertex has been used 
in a convex polygon * / 


if(track ptr->back->index == stack; TOP].index) 
track ptr->segment_ used= TRUE; 


/* for each vertex of the convex polygon, pass its 
coordinates to IRIS hardware with the point draw 
command */ 


count= I: 
while(count <= TOP) { 


while(track ptr->index != stackicount).index) { 
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track ptr= track ptr->forward: 


} 


if(track ptr->back->index == stack count- 1\.index) 
track ptr->segment used= TRUE: 

pdr(pts stack count .index) 0. 
pts'stackicount).index] 1 . 
pts'stack count .index| 2)): 

count+= lI: 


} 


/* pass the command to the IRIS hardware that the 
fill the completed convex polvgon *' 


pclos(): 
clear _stack(): 


return: 
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KOKO OO SO RIE OR IO KOK ak aa ak ok oe 


ROUTINE: obj complet.c 

BOO ORR Ra Ook kk dk dk kok ae 
This routine checks each vertex and edge of the original 
concave polygon to see if they have all been utilized in the 


process of decomposition. If so a value of | is returned. 
dO Ook kok kk kok dk OK OK kok kk a ok kkk kK ok aR ek ok 


int object complete(total pts) 


int total pts: 


pointe pls 
int count, value: 


count= 0: 
value= 1: 
pl= start ptr: 
‘* search through the linked list of structures. for the 
vertices, and return the index of the first vertex that 
creates an edge of the concave polygon not vet utilized * 
do { 
if(pl-> segment used) { 
pl= pl-> forward: 
count += 1]; } 


else { 
value= 0: 
} 
} while(count < total pts && value == 1); 


return (value): 
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DOK OOO Kk ka kk dak ok ok eo kee 


ROUTINE: intersecting.c 
OOK OOK EOE KOK a 
This routine takes two vertices and the slope and intercept 
of the line segment they create to determine if it intersects 
any of the edges of the original polygon. If so then a value 


of 1 is returned. 
OOK EO ROK RE EOS IOI OOO KK OK KK kkk ok 


intersecting lines(x2, y2. xl. yl. slope. intercept. total pts. pts) 


Coord x2. v2. xl, yi: 
float slope: 

Coord intercept: 

int total pts: 
Coord pts{|[3}; 


point “pl: 
float x intersect. y intersect. xX min. XxX max. y min, y max, fudge: 
int value: 


'* define a value to offset round off error * 
fudge= 0.000001; 


value= 0; 
pil = start. ptr. 


-* for each edge of the concave polygon we check if either vertex 
lies on the edge. If not then we determine the intersection 
of the line. defined by the edge, with the line defined by the 
two vertices. Finally we determine if this intersection lies on 
both line segments. If so then we have found an actual intersection. * / 


do { 
/* check if first vertex lies on the polygon’s edge being examined * / 


if(((x1 != pts|p1->index]{0)) | 
((x1 != pts|p1->back->index| 


(yl != pts|/p1->index]|1])) && 
[0}) || (y1 != pts[pl->back->index|[1]))) { 
/* check if second vertex lies on the edge */ 


if(((x2 != pts|pl->index![0)) || (v2 != pts|pl->index]|[1])) && 
((x2 != pts|pl->back->index!|0)) |; (y2 != pts|pl->back->index||1'))) { 
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/* if the lines are not parallel then ... * 
if(slope '= pl->slope) { 
/* if the line created by the vertices is vertical * 
if(slope == 99999.0) { 
/* determine the coordinates of the intersection * 
xX Intersect= x2: 
y intersect= pl->slope * x intersect+ pl->intercept; } 
/* if the edge of the polygon is vertical * / 
else if(pl->slope == 99999.0) { 
/* determine the coordinates of the intersection * 
x intersect= pts'pl->index| 0: 
y intersect= slope * x_intersect— intercept: } 


/* if neither line is vertical *— 


else { 
/* determine the coordinates of the intersection * 


x intersect= (intercept - (pl->intercept))/((p1->slope) - slope): 
y_intersect= slope * x intersect~ intercept: 


} 


/* determine the x range that the intersection must occur in 
to lie on the line segment created by the two vertices * / 


max min(x2.x1.&x max,&x_ min); 
/* check if the intersection falls within this range */ 
if((x intersect >= x_min-fudge)&&(x intersect <= x max+fudge)) { 


/* since the intersection is within this x range we now determine 
the x range for the edge of the concave polygon */ 


max min(ptsipl->index|/0], pts|pl->back->index||/0). &x max. &x_ min); 
/* check if the intersection falls within this range */ 
if((x_ intersect >= x_min-fudge)&&(x intersect <= x max-+fudge)) { 


/* now that the intersection falls within the x range of both 


line segments, we must calculate the y range for each 
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segment and perform the checks * 


/* y range for the vertices line segment * 


max min(y2,yl.&y max,&y min): 
if((y intersect >= y min-fudge)&&(y intersect <= y max~+fudge)) { 


/* vy range for the concave polygon edge *. 


max min(pts pl->index!:1 .pts pl->back->index!/1).&:y max.& y_ min): 
if((y_ intersect >= y_ min-fudge)&&(y intersect <= y max-+fudge)) { 
value= 1; 
/* the intersection exists for the actual line segments and 
therefore the vertex on top of the stack must be removed *, 


pl= pl->forward; 
} while((pl '= start ptr)&&(value == 0)); 


return(value): 


74 


4 OK ok OK kk OR dk bk dk bok ko okok kok dakook dock ko 


ROUTINE. max min.c 


KOK EO oR ok ek EO OR ROKK RIK QO OOK OKO OK EK ERO RO SOK x 


Given two float values, pl and p2. this routine determines 


which is larger of the two. 
ne OK KOR ROK ROR EE ROR ER OR ERROR OO EK OE ERS EE OE RK x 


max min(pl, p2. max. min) 


Coord pl. p2. *max,. * min: 


Mp) < p2) { 


*max= p2: 

*min= pl; } 
else { 

*max= pl: 

*min= p2: 


} 


return; 
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